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ABSTRACT 

In this follow-up paper, we continue our study of the effect of using knowledge from electromagnetic observations in the gravitational 
wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in 
the low-frequency range, 10~ 4 Hz < / < 1 Hz. In the first paper, we have shown that the strong correlation between amplitude and 
inclination can be used for mildly inclined binaries to improve the uncertainty in amplitude, and that this correlation depends on the 
inclination of the system. In this paper we investigate the overall effect of the other orientation parameters, namely the sky position 
and the polarisation angle. We find that after the inclination, the ecliptic latitude of the source has the strongest effect in determining 
the GW parameter uncertainties. We ascertain that the strong correlation we found previously, only depends on the inclination of 
the source and not on the other orientation parameters. We find that knowing the sky position of the source from electromagnetic 
data can reduce the GW parameter uncertainty up to a factor of ~ 2, depending on the inclination and the ecliptic latitude of the 
system. Knowing the sky position and inclination can reduce the uncertainty in amplitude by a factor larger than 40. We also find that 
unphysical errors in the inclinations, which we found when using the Fisher matrix, can affect the corresponding uncertainties in the 
amplitudes, which need to be corrected. 
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1. Introduction 

Even though population synthesis studies have predicted ~ 10 8 
galactic white-dwarf binaries (e.g. Nelem ans et al.|2001| l, there 
are only about ~ 50 compact binary sources that have been stud- 
ied at optical, UV, and X-ray wavelengths (e.g. Roelofs et al. 



2010). A space-based gravitational wave (GW) detector, such as 
the proposed ESA mission eLISA, is expected to observe mil- 
lions of these predicted compact Galactic binaries with periods 
shorter than about a few hours (Nelemans 2009, Ama ro-Seoane| 
|et al. |2TJT2| ), amongst other astrophysical sources, and resolve 
several thousands of these binaries (Nissa nke et al.| [2 012) . In a 
recent paper, we investigated the effect of using knowledge from 
electromagnetic (EM) observations in the gravitational wave 
data analysis of Galactic binaries by exploring the correlations 
that might exist between these parameters. We found that for 
binaries with relatively low inclinations (i < 45°), an EM con- 
straint on inclination can improve the estimate of the GW ampli- 
tude by a significant factor ( |Shah et al.|2012| Paper I, hereafter). 
The improvement depends on the signal-to-noise ratio (S/N) of 
the source and on the EM constraint on inclination itself. We 
also found that this correlation is a strong function of the incli- 
nation of the system, but that it is independent of other angular 
parameters that describe the binary, in particular, sky position, 
and polarisation angle. From the initial results we also found that 
several other correlations that were rather small as a function of 
inclination (e.g. Figure 4 in Paper I) change drastically when 
computed as a function of sky position. We address this problem 



here and also investigate the general dependence of parameter 
uncertainties on the orientations of the binary. Additionally, we 
discuss the influence of knowing the sky position from EM data 
on the parameter uncertainties of the binary sources. 

The angular resolution of cLIS Ap]for mo nochromatic bina- 
ries has been studied previously (Peterseim et al.|| 1997 Cutler 



T998l|Rogan & Bose|2006l|Biaut|2011| i. |Peterseim et al.| ( |1997) 
have shown the uncertainty in ecliptic latitude and ecliptic lon- 
gitude as a function of both the parameters for a specific case of 
a cross-polarised source radiating at / = 3mHz, which would be 
observable by cLISA. For a source with S/N = 115, depending 
on the values of the ecliptic latitude and longitude, they found 
the uncertainties in these angles to range from 1 to 8 millira- 
dians. However, they assumed that except for the sky position, 
all other parameters are known a priori. Their study was also 
limited by a very simplified detector response and was tailored 
for a system at an inclination which is favourable in extract- 
ing its parameters. Relaxing the assumption of prior knowledge 
of other parameters, [Cutler] (1998) has shown that for a binary 
for a few sets of sky positions, the solid angle uncertainty (ctq) 
ranges from 3 to 300 square degrees, depending on the frequency 
of the source. In subsequent work by Rogan & Bose (2006), 
sky position uncertainties were estimated for signals modelled 
for cLISA. They included chirp as an additional parameter, and 
these signals were not limited to a certain polarisation. Recently, 



1 classic LISA refers to the detector configuration in consideration 
prior to eLISA (see Paper I). 
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Blaut] ( |201 1) has shown sky position uncertainties for various 
combinations of data channels from cLISA and included de- 
tailed detector response for the two distinct cases of low- and 
high-frequency sources. For a 10 3 Hz source with an edge-on 
orientation and an S/N of 10, an observation of one year is pre- 
dicted to give a solid angle resolution of ~ 0.1 steradians which 
is a result very similar to the study by |Rogan & Bose] (|2006). 
However, none of these previous authors provided an overall ef- 
fect of varying multiple angular parameters, which can have a 
significant effect on parameter uncertainties. In this paper we 
aim to generalise the uncertainties in all parameters for all pos- 
sible orientations of the binary and calculate the improvement in 
these uncertainties that is possible thanks to prior knowledge of 
the sky position and/or inclination of the system. We show this 
for eLISA. 

A monochromatic binary expected to be observed by eLISA 
is completely described by seven parameters: dimensionless am- 
plitude (Jl), frequency (/), polarisation angle (if/), initial GW 
phase ((/>[)), inclination (cos t), ecliptic latitude (sin (3), and eclip- 
tic longitude (A). The inclination of the binary is defined as 
cos i = n-L, where n is the line-of-sight vector from the observer 
to the source and can be described by the sky position parame- 
ters f3 and A. L is the orientation of the orbital plane of the binary 
(or angular momentum vector), which can be specified by 6i and 
<p L . One of the main findings of Paper I is that the parameter un- 
certainties of a monochromatic binary and the three strongest 
correlations between these parameters (the correlation between 
amplitude and inclination cosl , the correlation between phase 
and polarisation c^, and the correlation between phase and fre- 
quency Off) are strong functions of the source's inclination t. 
Furthermore, we found that the parameter uncertainties of such 
a binary also change as a function of other orientation parame- 
ters. These are parameters that describe the relative geometry of 
the binary in relation to the detector: if/, A, and /3. It is natural 
to then ask the questions, how the possible strong correlations 
between the parameters due to these four orientation parame- 
ters influence these uncertainties, and to which extent knowl- 
edge from the EM data can improve the remaining parameter 
determinations. Of the remaining three parameters, J{ (see Eq. 3 
in Paper I) is intrinsic to the source and does not depend on its 
relative orientation to the detector. Since it only scales the S/N 
of the system, varying the amplitude thus will affect parameter 
uncertainties as expected, where larger J{ will give smaller un- 
certainties and vice versa. The initial phase </>o, which does not 
evolve for a monochromatic binary in the eLISA band, has no ef- 
fect on the signal. On the other hand, the effect of frequency on 
the parameter accuracies can be significant; the strength of the 
modulations in the signal (which is caused by the motion of the 
detector around the Sun and relative to the source) affects low-/ 
sources differently than high-/ sources. This has been addressed 
in detail in previous studies, e.g. Blaut] pOTT) . Hence, here we 
focus on the effects on the parameter uncertainties of a Galactic 
binary due to sky position and the orientation parameters. 

In this study, we consider low-frequency binaries. We start 
by briefly summarising the analysis methods and then describe 
our Monte Carlo (MC) simulations in Section 2. In Section 3, 
we summarise the important results and our interpretation. In 
Section 4 we show that uncertainties can be severely overesti- 
mated for parameters with clear physical bounds, such as incli- 
nation, and thus also for parameters that are tightly correlated 
with them. We also briefly discuss the S/N limit for which the 
results are valid in Section 4, and present our conclusions in 
Section 5. 



2. Monte Carlo simulations 

Our signal and noise modelling and our data analysis method 
have been described in Paper I and we briefly summarise them 
here. The signals from the binaries that are expected to be ob- 
served by eLISA are modelled as monochromatic sources char- 
acterised by seven parameters, as mentioned above. Most of 
these sources radiate at low frequencies (milliHertz), including 
some of the verification binaries like AM CVn, and at these low 
/, the detector response can be derived in the long-wavelength 
approximation. The noise is composed of instrumental noise de- 
scribed as a Gaussian random process, and a foreground noise 
mostly arising from the millions of double-detached Galactic 
white-dwarf binaries. After removing the binaries, which are 
predicted to be resolved due to their strong signals, the fore- 
ground noise is also described by a Gaussian distribution. Given 
that a signal has a known waveform and the noise is Gaussian, 
we can use the Fisher matrix to estimate the uncertainties in the 
parameters that describe the waveform. The inverse of the Fisher 
information matrix (FIM), known as the variance-covariance 
matrix (VCM), gives the uncertainties and the correlations be- 
tween the parameters. Given a number of parameters n, one can 
calculate the nxn VCM. For details we refer to Paper I. 

Uncertainties in parameters are proportional to the inverse of 
the S/N, but they are also influenced by the correlations between 
parameters. Generally, a system with higher S/N has smaller pa- 
rameter uncertainties and vice versa. However, if there is a strong 
correlation between any two parameters, this (in addition to the 
S/N) will influence their uncertainties. As mentioned in the pre- 
vious section, we find that the overall effect of the sky position 
and the orientation parameters of a binary has not been investi- 
gated and this is one of the goals of this study. To do this, we 
fixed the values of 3i, f, and tpo to those of AM CVn (see Paper 
I) and varied its if/, t, /3, and A within their allowed ranges us- 
ing a Monte Carlo method. We carried out a number of of MC 
simulations. Each simulation corresponds to varying a different 
(sub)set of these four parameters simultaneously. The motivation 
is to understand how each of the four parameters above influ- 
ence the GW data analysis individually and in combination with 
each other. We explain these motivations briefly below, which 
will become clearer in section 4 where the choice of each MC 
simulation is explained in more detail. 

First of all, we set the distance of the AM CVn system close 
enough so that the S/Ns are sufficient to believe the Fisher matrix 
results (see section 4). The first set of simulations we performed 
for an AM CVn of a fixed M, f, and <f>o, while its position and 
orientation parameters A, /3, if/, and i were all varied simultane- 
ously 4000 times. Thus, from the sample of these 4000 systems 
with various orientation parameters we first ascertained that the 
dominant effect on the parameter uncertainties is indeed the in- 
clination of the system (result from Paper I). Next, we can quan- 
tify how the other three angles affect the parameter uncertainties 
globally. We performed another set of calculations where we 
fixed inclination (to a few choices explained later), frequency, 
and the phase (like above), and we varied the amplitude (by 
changing the distance, d), the sky position, and the polarisation 
angle. For each choice of the inclination, MC simulations were 
performed for 500 systems to quantify the effect of varying the 
S/N (by changing the distance) and the ecliptic latitude, longi- 
tude, and the polarisation angle on the parameter uncertainties. 
Finally, to quantify the possible improvement in the parameter 
uncertainties gained by prior (EM) knowledge of the sky loca- 
tion, we carried out 500 MC simulations for an AM CVn at a 
fixed inclination and varied the rest of the orientation parameters 
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and the distance. However, instead of a full 7x7 VCM, we calcu- 
lated a 5 X 5 VCM, where for each system the randomly chosen 
sky locations were fixed to that value. Similarly, we repeated this 
exercise to quantify the improvement in the parameter accuracies 
gained by EM knowledge of inclination in addition to sky posi- 
tion, and thus we calculated a 4 x 4 VCM for each system. For 
clarity, we briefly summarise these four sets of MC simulations 
below: 

1 . MC 1 : Fix J[ (i.e. d,m\, and m.2 where, w, is the mass), /, and 
0o- Randomly pick the orientation of the binary, L{(pL,0L), 
A, and /? over their physical ranges. The parameters if/, and t 
are then given by the following equations ( |Cornish &"R ubbo 
|2003) : 

cost = cos(#£) sin(/8)+sin(#£) cos(/?)+cos(0£-/l) and (1) 

sin/? cos(/t - <b£) - cos(# L ) sin(6» L ) 

tan ib = . (2) 

sin(0 L ) sm{A-<f> L ) 

2. MC2: For each i = 30°, 60°, 90°, randomly pick A, f3, if/, and 
d. 

3. MC3: Like MC2, but fix A, and /? assuming an a priori known 
sky position [5x5 FIM]. 

4. MC4: Like MC2, but fix A, f3, and ( assuming an a priori 
known sky position and inclination [4x4 FIM]. 

Each of the four parameters are randomly picked from their uni- 
form distribution in the following range: A e [0,2^], sin/? e 
[-1,1], <f) L e [0,2n], 6 L e [0,n], and d e [0.01,0.3] kpc. 
The corresponding ranges in t and if/ are [0, 7r]. We also define 
the following inclinations for presenting the results below (un- 
less specified): face-on:= 0° < i < 45° or 135° < t < 180°, 
mildly face-on:= 45° < i < 60° or 120° < l < 135°, mildly 
edge-on:= 60° < i < 80° or 100° < l < 120°, and edge-on:= 
80° < i < 100°. 
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Fig. 1. Two-dimensional distributions of signal-to-noise ratio as 
a function of cost, sin/?, A, and if/. Face-on, mildly face-on, 
mildly edge-on, and edge-on systems (whose exact bounds are 
specified in Section 2) can be read off from the top-left panel. 
In the other panels, the trend from face-on to edge-on systems 
follows from the top of the distribution to the bottom. 



3. Results 

3.1. Global dependence of the uncertainties: dominance of 
inclination (MC1) 

We have seen in Paper I that the parameter uncertainties are a 
strong function of the inclination of the system, which is due to 
the strong dependence of the signal strength, i.e. S/N, on the in- 
clination and to the strong correlations. In this section, we first 
discuss the global behaviour of the parameter uncertainties as a 
result of varying the sky position, inclination, and the polarisa- 
tion angle at the same time. We verify the result from Paper I that 
inclination has a dominant effect on the parameter uncertainties. 
We show that after the inclination, the ecliptic latitude, /?, of the 
system has the strongest influence on the parameter uncertain- 
ties. 

Using the results from MCI, we show in Figure [T] the S/N 
of binary systems as a function of the four angular parameters 
over which they are randomly picked. As expected, there is a 
strong dependence of S/N of a system on inclination, as shown 
in the top-left panel. Due to this dependence on the inclination, 
in the other panels we find that the face-on systems form a band 
at the top of these plots, followed by bands of mildly face-on, 
mildly edge-on, and edge-on systems, in that order. In the top- 
right panel, there is a mild but noticeable increase in S/N, for 
the binaries at the ecliptic plane, sin/? = 0, whereas the eclip- 
tic longitude and the polarisation angle have no significant effect 
on the binaries' S/N, as shown in the bottom panels. Figure [2] 
is a plot of parameter uncertainties as a function of sin /? taken 
from MCI. In the top panels, uncertainties in M, i, <f>, and if/ are 
smallest for edge-on systems (forming a band at the bottom) and 
progressively increase for the face-on systems (the top of distri- 
bution). This is discussed in Paper I and results from the corre- 
lations between these parameters. The lower inclination systems 
t < 45° or t > 135° have very similar signal shapes in those 
ranges, whereas systems with close to edge-on orientations are 
distinguishable by both the shape and structure for small dif- 
ferences in inclinations. For face-on systems a small change in 
inclination is indistinguishable from a small change in its intrin- 
sic amplitude, whereas for an edge-on system, a small change 
in inclination produces a noticeably different signal. Thus even 
for their lower S/Ns, the inclinations of the edge-on systems are 
better determined than those of the face-on systems. A similar 
argument based on the strong correlation between <p and if/ for 
face-on systems can explain the huge uncertainties of these pa- 
rameters despite their strong signals (see Paper I). Since there 
are no such correlations between /, A, and /? (or the solid angle 
Q), their uncertainties should be a consequence of the S/N, i.e. 
from the top panels of Figure [T[ a higher S/N implies lower val- 
ues for the uncertainties and vice versa, which indeed holds for 
(Tf and a a- From Figures [T|and[2]it is evident that inclination has 
the strongest effect on the S/N and so to the leading order this pa- 
rameter affects the parameter uncertainties, which confirms our 
results in Paper I. 

After inclination, the strongest effect on the S/N of the sys- 
tems is produced by the ecliptic latitude. Thus, we expect the 
uncertainties in the parameters also to be affected by /? (see top- 
right panel in Figure [TJ; sources at the ecliptic poles should have 
larger uncertainties than those at the ecliptic plane. In Figure [2] 
we can see that for a fixed distance to the source d, most un- 
certainties {ctji, cr t , O-0, cr^j, cry, and cr^) are indeed smaller to- 
wards the plane than close to the poles. This is probably largely 
due to the S/N as a function of sin /?. The unusual behaviour of 
cr^, which increases at the plane compared to its value at the 
poles (sin/? = ±1), has been pointed out in |Cornish & Larson] 
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Fig. 2. Two-dimensional distributions of the uncertainties in the seven GW parameters of the galactic binary systems as a function 
of sin/3 from MCI as described in the text. Units of the angles are in radians. 
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Fig. 3. Uncertainties in all seven parameters of the systems as a function of S/N from MC2. Open circles, squares, and filled triangles 
represent binary systems at fixed inclinations of 90°, 60°, and 30° respectively. The size of the symbol represents the magnitude of 
\/3\ for each simulated binary. In the bottom panels the different groups of symbols overlap. 



(2003} and is caused by the Doppler modulation Oq (Eq. 9 in 
Paper I). This behaviour is surprising since O d has the largest 
magnitude on the plane (since it is proportional to cos f3), and so 
we naively expect the sky resolution to be better on the ecliptic 
plane. However, the resolution in j3 is determined by the rate of 
change of the Doppler modulation with respect to /3 (i.e. propor- 
tional to sin p), which is the highest at the pole and lowest at the 
plane. This is also evident in the solid angle uncertainty cr n in the 
bottom-right panel in Figure[2] where ctq - 27r(cr sin/3 cr /1 -Csin J s,i), 
Csingi b eing the unnormalised correlation between sin/3 and A 
( |Cutler|[T998[ l. We also observe that cr^ has a distinct U shape, 
becoming very severe at the poles. This is merely because A be- 
comes progressively less well-defined as a source lies closer to 



the pole. Note that at the pole A is not defined, which means that 
it can take any value. 



3.2. Dependence of the uncertainties on ft: S/N or 
correlations (MC2)? 

The general behaviour of the parameter uncertainties (except erg, 
hence cr n ) in Figure[2]shows a minimum at the ecliptic plane that 
could be attributed to the S/N, since it is the highest for sources 
in the plane. To investigate this change in the uncertainties as 
a function of sin/3, we consider the results from MC2 and the 
correlations between parameters from MC 1 . In MC2 we varied 
S/N (by varying the distance to the binary, d), A, {!, and \p for a 
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Fig. 4. Two-dimensional distributions of the normalised correlations as a function of sinyS from MC2. We show here only the strong 
(normalised) correlations, meaning |cy| > 0.5. 



few fixed inclinations. We fixed the t for each of the MC2 simu- 
lations because we established its dominant effect through the 
S/N in the previous section, therefore the goal now is to de- 
termine which contributes most to the behaviour observed in 
Figure [2] the S/N (top-right panel of Figure [TJ, or the correla- 
tions as a function of sin If the larger uncertainties at the pole 
compared to the plane are solely due to the corresponding S/Ns, 
then we should see from the MC2 simulations that uncertainties 
anti-correlate tightly with the S/N of the sources. In Figure[3]we 
show three distinct curves in S/N vs. parameter uncertainty that 
roughly show monotonic functions where all the uncertainties 
decrease with increasing S/N. However, all uncertainties except 
cry have a significant horizontal scatter. In the figure, the size 
of the symbol represents the magnitude of the ecliptic latitude. 
Higher values of \/3\, which represent regions around both ecliptic 
poles, concentrate towards the right-hand side of the plots (with 
the exception of erg), and the opposite is true for the lower values 
of 1/3*1, which roughly represent the ecliptic plane. In the top pan- 
els, from left to right, the curves correspond to l = 90°, 60°, and 
30° respectively, and these overlap in the bottom panels. As ex- 
plained above, cr^, cr L , cr$, and cr^ are larger for face-on systems 
(i = 0,30°) than for edge-on systems [i = 90°), owing to the 
correlations between these parameters even though face-on sys- 
tems have higher S/Ns. Since these correlations are very weak 
for A, p, and i/r for all inclinations (see Paper I), their uncertain- 
ties overlap in the bottom panels of Figure[3] This spread across 
the horizontal axis must be due to changes in their correlations 
as a function of the ecliptic latitude. This is indeed shown in 
Figure |4j which is taken from MC 1 . In the figure we show a 2d 
histogram to display the global dependency of the correlations 
on the ecliptic latitude. We notice that the correlations involving 
A, i, <p, and J{ (cyu> c <pn C <M> c iA, and c^) are generally weaker 
at the ecliptic plane. However, at the poles these correlations can 
take any value from to ±1 and concentrate at +1. This is the 
case for all inclinations. Accordingly, close to the ecliptic poles, 
there are sets of A and iff for which the correlations can become 
very strong, which explains the horizontal spread of uncertain- 
ties in Figure [2] Note that in Figure [4] we have shown (only) the 
normalised correlations cu, whose absolute values can exceed 



0.5. A detailed explanation of these correlations can be found in 
the Appendix. 

3.3. Gain from a priori knowledge of sky position and 
inclination (MC and MC4) 

From Figures [2] and [3] we observe that uncertainties in sky loca- 
tion are rather large (which has also been pointed out in previous 
studies). For the S/Ns in this study (~ 8 — 100), the uncertain- 
ties in crn range from ~ 0.2 to 6 square degree, as shown in 
Figure [Tj If EM data can provide a more accurate sky position 
and possibly an inclination for AM CVn system, we can use 
these independent (EM) measurements to improve accuracies in 
other parameters. To see by which factors the uncertainties im- 
prove, we fixed the sky position of each source and calculated 
the corresponding 5x5 VCM. The ratios of uncertainties from 
using a 7 x 7 VCM to using a 5 x 5 VCM, i.e. the improvement 
when knowing the sky position (ecliptic latitude and longitude) 
exactly, are shown in Figure [5] calculated from MC2 and MC3. 
In the figure, the size of the symbol represents the S/N of that 
system. The figure shows that the improvement factor depends 
on the inclination of the source: for edge-on systems (in open 
circles with t = 90°) cr^ can improve by a factor of ~ 3, but for 
systems closer to a face-on configuration in squares and filled 
triangles the improvement can only be up to a factor of ~ 2.2. 
One useful result from this figure is that cr^ and cr, can improve 
by a factor of ~ 2 for the systems close to the ecliptic poles. This 
gives a motivation for independent EM observations of the sky 
location. However, for edge-on systems there is no improvement 
in <Tj{ even though cr L improve^] This is simply due to the lack 
of correlation between these two parameters for l = 90°. 

In Paper I we found that for the HM Cnc system, the im- 
provement of cr^ by using the EM constraint in ( was smaller 
than for AM CVn, although the two have similar inclinations 
(thus the same (normalised) correlation between J[ and l). This 
was due to the relatively high S/N of HM Cnc compared to 
the latter and thus the individual improvement factors shown in 



2 Depending on the precision of the GW accuracy of the inclination, 
this may aid in finding eclipsing binaries. 



5 



Shah et. al.: Improving GW parameter estimation for compact binaries with LISA 







l.lOr 




1.08- 


■< 


1.06- 


< 


1.04- 








1.02- 












1.00- 




0.98 - 




0.96 - 



-1.0 -0.5 


0.0 0.5 1.0 


□ 




□ 
ft 


□ o 






- ■& 

















0.0 

sin/3 



■1.0 



-0.5 



0.0 

sin/3 



0.5 



1.0 



Fig. 5. Uncertainty ratios of the parameters between analyses using a 7 x 7 VCM and a 5 x 5 VCM, where the sky position is fixed in 
the latter, as a function of sin/3 (from the MC2 and MC3 simulations). Filled triangles, squares and open circles represent systems 
at fixed inclinations of 30°, 60° and 90° respectively. The size of the symbol for each binary represents the S/N of the system. 
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Fig. 6. Same as in Figure Hbut now we assume that in addition to sky position we also know the inclination of the binary. The 
uncertainty ratios of the parameters are between 7x7 VCM (MC2) and 4x4 VCM (MC4). The improvement factors in amplitude 
include corrections due to overestimated uncertainties in the amplitudes of face-on systems (see Section 4. 1) Note that the last panel 
is plotted in log scale in order to highlight the relatively large improvements for the systems with smaller inclinations. 



Figures [5] and [6] will also depend on the individual S/N of the 
systems (see Paper I). 

If in addition to sky position the inclination can be con- 
strained, as discussed previously in Paper I, the improvement 
factor in amplitude can be much larger because of the decou- 
pling of the correlation between amplitude and inclination. We 



show the ratio of parameter uncertainties from calculating the 
7x7 VCM to the uncertainties from a 4 x 4 VCM (where A, (3 
and t are fixed) in Figure [6] The improvements in <x^ and cr^ are 
roughly similar to the case in Figure [5] A remarkable difference 
is seen in the improvement of cr^ for binaries with l = 30° rep- 
resented in filled triangles plotted in the bottom-right panel of 
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Figure [6] This is expected because by fixing l, the strong cor- 
relation between amplitude and inclination is decoupled. Hence 
for the higher S/Ns of close to face-on systems, the uncertainties 
in 3i is significantly reduce. Note that here we assumed that the 
inclination is known exactly. 

4. Discussion 

4.1. Limiting inclinations to their physical range 

Some of the uncertainties in l in Figure[2]from MCI are notewor- 
thy in that they exceed their physical bounds, for example <x, =7 
radians. The cause for this is that the linearised signal approxi- 
mation that is internally employed in the Fisher matrix studies is 
not sensitive to the boundecFlparameters that describe the model 
signal, as has been pointed out in Vallisneri (2008). Thus, when 
an uncertainty in a (cyclicly) bounded parameter like the inclina- 
tion exceeds its physically allowed range, it means the quantity 
cannot be determined from GW data analysis alone. For the pa- 
rameter i, we can further ask what the unphysical uncertainty in 
l implies for the uncertainty in J[ for the mildly inclined systems 
that has a strong correlation between the two parameters. Since 
the Fisher matrix estimates cr, to unphysical values while in real- 
ity it is constrained within [0°, 90°], does this mean that the un- 
certainty in the amplitude (due to the correlation with i) is also 
overestimated? Or should we not trust the uncertainty in all those 
parameters where the cr, becomes unphysical? Consider the M-i 
error ellipse for a given signal h, with true signal values to and 
J?lo, where the two parameters are strongly correlated and where 
cr, >> n. The corresponding accuracy cr^ is then also highly in- 
flated. Now consider a second signal h{&li,i{), inside the error 
ellipse, but far from the true values, and such that L\ — i$ + 2nn 
with n > 1. From its cyclic behaviour, h(jM\,i\) = 
should hold. However, this last point would typically lie far out- 
side the error ellipse in the 3i - i plane, which would indicate 
that this is actually a signal that is very dissimilar to the true 
signal. Hence, the fact that the inclination is physically limited 
to the range [0, n] and that this is not taken into account in the 
FIM method means that we should not only limit the uncertainty 
in inclination to its physically allowed range, but that we should 
limit the uncertainty in amplitude accordingly. Specifically, the 
corrected amplitude uncertainty, cr^, is given by, 



o-' m = cr^[VCM] x 



cr t [VCM]' 



(3) 



where cr,' is the corrected inclination uncertainty. 

The value of inclination where its uncertainty becomes un- 
physical can occur at two places, i.e. where i + cr, > and/or 
l - cr, < 0, as shown in the left panel of Figure [7] In the panel, 
the 1-cr uncertainty ranges of the inclination are plotted against 
the sources' inclination for two cases of constant S/N of 10 and 
53. The thick horizontal dashed and dashed-dotted lines indicate 
the conditions where cr, becomes unphysical. At a fixed S/N of 
10 (in plusses), these values are at l ~ 48° and l ~ 47°, respec- 
tively, as indicated by the open circles, whereas for a higher S/N 



3 Note that choosing cos ( instead of t does not alleviate this problem 
since cos ( is also a bounded quantity and there is no formalism in FIM 
to handle bounded quantities, whether cyclic or not. In fact, for the FIM 
calculations we evaluate uncertainties in cos i and sin/?, while in quoting 
our results they are converted to uncertainties in i and /? by scaling them 
as cr, = cr cos ,/ s'mi and erg = cr sin ^/ cos/?, since these quantities are 
more intuitive. 

4 Since t is symmetric about n/2, we only considered the inclinations 
in the range [0, n/2]. 



of 53 (in stars) the above conditions for cr, occur at l ~ 27° and 
t ~ 32° indicated by filled circles. We calculated the inclinations 
at which these intersections occur for a range of S/Ns, 10 - 1 10 
and show the results in the right panel of Figure [7] These curves 
can be used to apply corrections to the limits of cr, in Figure [2] 
from MCI, which in turn can be used to limit cr^ from Eq. 3. In 
Figure|8]we compare the uncertainties in J[ corrected in this way 
with those from Figure [2] The corrected values for cr^ differ for 
the high S/N systems, as can been seen by comparing the two 
right-hand panels in Figure [8] In particular, the tail at the large 
crnIJi present in the top panel is absent from the bottom panel. 

There is a caveat to applying these corrections. In the left 
panel of Figure [JJ we can understand the cut-off at 0° at the bot- 
tom of the figure. However, at the top, the cut-off of 90° is not so 
clear. A system with a lower inclination (e.g. 20°) cannot have 
an error larger than ~ 50° because a system with an inclination 
in the range ~ 80° - 90° is distinguishable from a system with 
an inclination of ~ - 60°, as has been pointed out in Paper 
I. Thus the cut-off of i + cr, = 7r/2 for lower inclinations is a 
conservative estimate. The estimated posterior probability distri- 
bution functions (PDFs) obtained using the FIM do not predict 
the real PDFs for lower inclinations and lower S/Ns. Thus a full 
Bayesian analysis that takes into account the limits of the cyclic 
parameters as priors is needed to fully determine the PDF. This 
is beyond the scope of this paper, and we used these (probably) 
over-estimated uncertainties to give the corresponding bounds in 
the amplitude. 
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Fig. 8. Corrections applied to the uncertainties in amplitude and 
inclination from Figure [2] The row titled 'new' was produced 
using the curves from the right panel of the Figure [JJ on MCI 
simulations. 



In the same manner that cr^ was corrected in Figure [8] we 
corrected for the MC2 calculations. The predicted improvement 
factors of the amplitude in the bottom-right panel of Figure [6] 
already include the corrections mentioned above. 

4.2. Use of FIM at S/N ~ 15 

Another problem regarding the trustworthiness of the Fisher ma- 
trix is the limit of the S/N at which one can trust the values from 
VCM matrix. This has been addressed in Nicho lson & Vecchiol 
( 1998 ), where the authors have performed Fisher matrix studies 
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Fig. 7. Left: 1-cr uncertainty in the inclinations vs. source inclinations for two values of S/N: plusses (outer curves) correspond to 
S/N=10, stars (inner curves) correspond to S/N=53. The curves above the dotted line correspond to i + cr,, and those below to l - <r L . 
The thick horizontal dashed and dashed-dotted lines are the physical limits of t: below the line of t ± cr, =0° and above the line of 
l ± <x, = 90°, cr, can be truncated, since that value of ( is bounded between [0°, 90°]. Right: 'Unphysical C versus S/N for the systems 
from MCI. The thick dashed-dotted line calculated from the condition l - cr l = 0°, and the thick dashed line from the condition 
l + cr, = 90° demonstrated in the left panel. At an S/N of 10, the in open circles show a lower limit and the filled circles show an 
upper limit at the intersections with the two curves from the conditions above. 



in predicting the uncertainties of the parameters of the coalesc- 
ing binaries. They have compared these predictions with calcula- 
tions from Bayesian uncertainties, and their main finding is that 
above an S/N of 15, the uncertainties from the Bayesian calcu- 
lations converge to Fisher matrix predictions. At an S/N of 10, 
the analytic predictions of Fisher matrix deviate by about 6% 
and at a lower S/N of 7, they increase to about 25%. For MCI, 
the lowest S/N we considered is above ~ 15, as indicated in 
Figure [T] and consequently the results following MCI are rea- 
sonable. However, real observations for eLISA will of course 
also contain lower S/Ns. 



5. Conclusion 

We presented a follow-up study of correlations between param- 
eters that characterise Galactic binaries that will be observed by 
eLISA in the frequency range 10~ 4 Hz < / < 1 Hz. In Paper I 
we investigated the correlations between parameters by consid- 
ering three of the verification binaries (which are already known 
optically and are thus guaranteed sources for eLISA), and we ex- 
plained these correlations and parameter uncertainties as a func- 
tion of the inclination of the system. In this paper, we addressed 
the effect of all position and orientation parameters on the un- 
certainties and the possible strong correlations between them. 
We found that the overall effect is dominated by sky position 
and explained that dependence. Our main conclusions are: 

1. After inclination, the ecliptic latitude of the source has the 
strongest influence on its parameter uncertainties, mainly 
due to its dependence on the signal-to-noise ratio and the 
correlations that can become very strong at the poles. For 
a given source distance, the parameter uncertainties (except 
for erg) are larger for a source at the ecliptic pole than for a 
source at the ecliptic plane. This is because a signal received 
from a source at the plane is stronger than that from a source 
at the pole. 



2. In general, the uncertainty in sky location ctq is large, i.e. in 
the order of square degrees. By constraining A and /? from 
optical data, the uncertainties in ifr, </>, l, and J[ in the GW 
analysis can be improved by a factor of ~ 2. This factor is 
dependent on the inclination and the ecliptic latitude of the 
system. 

3. By constraining t in addition to the sky location, the improve- 
ment in <Tj{ is remarkably large (factors of 10 depending on 
the S/N of the source). Systems with t = 30° can have im- 
provement factors as high as 60. 

4. The remaining angular parameters A and iff do not influence 
the correlations or uncertainties significantly. 

5. The uncertainties in amplitude for low-inclination systems 
for which the corresponding uncertainties in the inclinations 
are unphysical (i.e. beyond their physical bounds) are also 
exaggerated. Fisher matrix predictions for these uncertain- 
ties have to be corrected for such systems. 
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Appendix A: Correlations as a function of p 

We found that several correlations between source parameters 
become amplified at certain ecliptic latitudes and thus are a 
strong function of /3, as shown in Figure |4] Several of these cor- 
relations do not have preference on the inclination of the system, 
and therefore we show a 2d histogram to show the dependency 
on the ecliptic latitude. In the figure we show (only) the corre- 
lations that become strong, i.e. where the normalised correlation 
|cy| > 0.5. Some of these correlations have been explained in 
Paper I and are only dependent on the inclination of the system; 
we show them here to confirm that claim. The remaining correla- 
tions and their dependency on the latitude can be understood by 
considering the orientation of the binary relative to the detector. 
We comment on or explain each of them below, 

1. cj\ L : Explained in Section 3.1. 

2. and c^f. Similarly to c^, these are strongly influenced by 
the source's inclination. The correlation in the bottom- 
left panel has a distribution that looks very similar to the case 
of Cft L and here too, phase and polarisation are strongly cor- 
related for face-on and mildy face-on systems, since these 
two angles describe similar properties. For edge-on systems, 
the two angles are defined on two planes that are perpen- 
dicular to each other, making them distinct. The correlation 
cqj in the top-right panel is strong for edge-on systems be- 
cause a small shift in phase is very similar to a small change 
in frequency. However, for face-on systems, the phase and 
frequency become degenerate owing to the uncertainty in 0, 
which is inflated due to the strong correlation with if/. There 
is a slight dependency of this correlation on sin/J for the 
edge-on orientations since the correlation has a larger spread 
at the ecliptic poles. This is due to the spread of the Cm for 
t > 80° systems where the spread in the phase at the poles 
also inflates the frequency. However, t < 60°, for has 
no such spread, and thus the distribution of vs. sin/? for 
these systems lies in a narrow band throughout. 

3. c^x- This can be understood by considering a binary with 
i = 90° close to one of the poles; rotating the binary about 
its line-of-sight is equivalent to changing its if/ and A because 
they are defined on the same plane (perpendicular to the line- 
of-sight). This also holds for a face-on binary, and thus if a 
source is close to the poles, if/ and A can be highly correlated. 
The spread of this correlation is most likely due to the fact 
that A is not well defined at the poles and thus the inflated 
uncertainty in the longitude implies a spread of correlations 
at the poles. 

4. c^y. This correlation is explained by the combination of c^x 
and c^, which are explained above. Since if/ and A can be 
highly correlated at the poles and <p and if/ are strongly cor- 
related for face-on systems, we found that <p and A should 
also be highly correlated at the poles, at least for the systems 
with (mildly) face-on orientations. For the edge-on systems, 
this strong correlation can also exist because even though 
C0,/, is not as strong in this case, the values of this correlation 
(as pointed out above) spread at the poles for i = 90°. This 
spread is responsible for the spread of cm around the poles 
for edge-on systems. 

5. ca and Csu- A strong correlation between inclination and the 
longitude is intuitively not as clear. However, this correlation 
combined with c^, for face-on systems would explain the 
correlation between amplitude and longitude, c^x- 

6. Cm. Like above, the correlation between phase and inclina- 
tion follows from the combination of c^x. an d c a, both of 
which are stronger close to the poles. 



